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Abstract 

In this paper, we investigate the performance of bidirectional relay selection using amplify-and- 
forward protocol with imperfect channel state information, i.e., delay effect and channel estimation error. 
The asymptotic expression of end-to-end SER in high SNR regime is derived in a closed form, which 
indicates that the delay effect causes the loss of both coding gain and diversity order, while the channel 
estimation error merely affects the coding gain. Finally, analytical results are verified by Monte-Carlo 
simulations. 

Index Terms 

bidirectional relay selection, analog network coding, imperfect channel state information 

I. Introduction 

Bidirectional relay communications, in which two sources exchange information through intermediate 
relays, have gained a lot of interest by now, and different transmission schemes have been proposed 
All . In (21, 0, an amplify-and-forward (AF) based network coding scheme, named as analog network 
coding (ANC), was introduced. With ANC, the data transmission process can be divided into two phases, 
and the spectral efficiency, which is restricted by half-duplex antennas, can get improved. Recently, relay 
selection (RS) for bidirectional relay networks has been intensively researched to achieve full spatial 
diversity and better system performance, which requires fewer orthogonal resources in comparison of 
all-participate relay approaches H, Q. Performing RS, the "best" relay is firstly selected before data 
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transmission by the predefined criterion floTl- lfTTIl . In 0, Q, the authors proposed the max-min sum rate 
selection criterion for AF bidirectional relay. In lf8Tl — iTTTTl . selection criterions in minimizing the symbol 
error rate (SER) were introduced and analyzed. 

To the authors' best knowledge, most works about RS in bidirectional relay only consider perfect 
channel state information (CSI). However, imperfect CSI, i.e., delay effect and channel estimation er- 
ror (CEE), has great impact on the performance of bidirectional relay selection. Specifically, the time 
delay between relay selection and data transmission causes that the selected relay may not be optimal for 
data transmission lfT2l - |[T4l . And similarly, channel estimation errors can not be ignored either lfl31 - lfl"8l . 
In |fT8l , the authors analyzed the performance loss of bidirectional relay selection using decode-and- 
forward protocol with CEE, but the impact of imperfect CSI on a general bidirectional AF relay selection 
was not provided. 

In light of the aforementioned researches, we analyze the impact of imperfect CSI, including delay 
effect and CEE, for bidirectional AF relay selection in this paper, which has not been studied previously. 
The asymptotic expression of end-to-end SER is derived in a closed form, and verified by computer 
simulations. Analytical and simulated results reveal that delay effect reduces both the diversity order and 
the coding gain, while channel estimation error merely causes the coding gain loss. The main contribution 
of this paper can be summarized as follows: 

1) The asymptotic SER expression for bidirectional relay selection is provided in a closed form, which 
matches the simulated results in high SNR regime; 

2) Imperfect CSI, i.e., delay and channel estimation errors, is taken into account to derive the analytical 
results, and its therein impact is investigated. 

The remainder of this paper is organized as follows: In Section II, the system model of bidirectional AF 
relay selection, and the imperfect CSI model are described in detail. Section III provides the analytical 
expression of bidirectional relay selection with imperfect CSI. Simulation results and performance analysis 
are presented in Section IV. Finally, section V concludes this paper. 

Notation: {■)* and |-| represent the conjugate and the absolute value, respectively. E is used for the 
expectation and Pr represents the probability. The probability density function and the cumulative 
probability function of variable x are denoted by f x (•) and F x (■), respectively. 

II. System Model 

The system investigated in this paper is a general bidirectional AF relay network with two sources 
Sj, j = 1,2 exchanging information through the intermediate N relays E4, i = 1,...,N. The direct 
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link between Si and S2 does not exist, and each node is equipped with a single half-duplex antenna. 
The transmit power of the sources is assumed to be the same, denoted by p s , and all the relays have 
the individual power constraint, denoted by p r . The channel coefficients between sources and relays are 
reciprocal, and these coefficients are constant over the duration of one data block. 

The whole procedure of bidirectional AF relay selection is divided into two parts periodically: relay 
selection process and data transmission process, which will be described concretely in the next section. 
Let h s ji and h S) ji represent the actual and the estimated channel coefficients between Sj and Ri during 
the relay selection process, respectively; let htji and htji represent the actual and the estimated channel 
coefficients between Sj and Ri during the data transmission process, respectively. All the actual channel 
coefficients are independent identically distributed (i.i.d.) Rayleigh flat-fading with zero mean and unit 
variance, i.e., E (\htji\ 2 ^ = E ^\h s ji\ 2 ^j = 1, and thus, \h s ji\ 2 and \h tt ji\ 2 are both exponentially 
distributed with unit mean. Both the sources can know the global channel coefficients by estimating the 
training symbols, while each relay only has its local channel information. 



A. Model of Delay Effect 

Due to the time delay between relay selection process and data transmission process, h s ji is not the 
same as htji, which means the CSI is outdated. Their relationship can be modeled by the first-order 
autoregressive model lfl4l : 



Kji = Pft Ka + v 1 " PfCj , (!) 

where htji is a zero mean complex-Gaussian RV with variance of o\ t ; h s ji and sj are i.i.d. random 
variable (RVs) with zero mean and variance of a 2 h and a\., respectively. In this paper, we assume 

The correlation coefficient pj ( < pf. < 1, where pj = 1 represents no delay effect, in other words, 
the CSI is not outdated) between Sj and relays is defined by Jakes' autocorrelation model |[T4ll : 

Ps . = J {2irf dj T) , (2) 

where Jo (•) stands for the zeroth order Bessel function E3l . fa is the Doppler frequency, and T is 
the time delay between the relay selection process and the data transmission process. In this paper, 
two variables pp, j = 1,2 are used to represent the correlation coefficients between Sj and the relays, 
respectively, for f dl and fa 2 may be different. 
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B. Model of Channel Estimation Error 

Let h denote the actual channel coefficient and h represent the estimated channel coefficient, and then 
their relationship can be modeled as follows lfT5l : 



and 



h = h + e , 



h = p e h + d 



(3) 



(4) 



where h and CEE e are independent complex-Gaussian RVs with zero mean and variances of a\, a 2 , 
respectively, h and CEE d are also independent complex-Gaussian RVs with zero mean and variances 
of <t?, a 2 D , respectively. The correlation coefficient p e = o~\l&\ (0 < p e < 1. where p e = 1 means 
no CEE) is determined by the concrete channel estimation method. In addition, p e can be modeled as 
an increasing function of the training symbols' power P, i.e., p e — > 1 when P approaches infinity |fT9l , 
EOl . In this paper, we assume p e = Pj (P + N ) ED. 

According to the above relationship, the variances of CEE are given by : 

.2 \ 2 



^ 2 J2 



1 , i '1 ^ Pe 2 
°h = I 1 - Pe) °~ h = — <?h 



Pe 



and 



l l 2 2 /i \ 2 

(7 D = (T h- Pe a h = I 1 - Pe)0- h 



'h ■ 



(5) 



(6) 



Assuming a\ = 1 in this paper, we have <r? = p e 1 and a 2 D = 1 — p e according to ((5]) and ©. 



C. Relationship between h s ji and htji 

For the bidirectional relay selection communications, h s ji is used for relay selection, and h t ji is used 
for data detection. According to the model of imperfect CSI, we have : 
Lemma 1: htji and h s ji can be related as : 



hji = Pjh s ji + J I - ppj , 



(7) 



where Vj and h s ji are i.i.d. RVs, and 



Pi 



i, if ph = !; 

PeP f f , if Pfj < 1. 



When the CSI is not outdated, i.e., pf = 1, j t ,ji 
distribution. 



H,ji 



and %ji 



"S,]l 



(8) 



have the same 
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When the CSI is outdated, i.e., pf. < 1, the probability density function (PDF) of ^ t ,ji conditioned by 
%ji can be expressed as : 



faM^i (y\ x ) 




(9) 



where Iq(-) stands for the zeroth order modified Bessel function of the first kind [23], and 

h h 3 ji 

Proof: The proof of Lemma 1 can be found in Appendix A. 



III. Performance Analysis of Bidirectional Relay Selection with Imperfect CSI 
A. Instantaneous Received SNR at the Sources 

As mentioned above, the whole procedure of bidirectional relay selection is divided into relay selection 
process and data transmission process. 

In the relay selection process, the central unit (CU), i.e., Si or 5*2, estimates all the channel coefficients 
h s ji. Then, based on the predefined selection criterion, CU selects the "best" relay from all the available 
relays for the subsequent data transmission and other relays keep idle until the next relay selection instant 
comes. There are several selection criterions for bidirectional relay lf6ll — iTTTTl - In this paper, we adopt the 
Best-Worse-Channel method for relay selection which has the best performance in minimizing the average 
SER and is tractable for analysis iPTOl , ifTTll . According to this criterion, the index k of the selected relay 
satisfies : 



k = arg max min 

i 



h s ,n 



s,2i 



and thus, 



mm 



s.lk 



max mm 



s.li 



s,2i 



(10) 



(11) 



The subsequent data transmission process can be divided into two phases. During the first phase, 
the sources simultaneously send their respective information to the intermediate relays where only the 
selected relay is active. The superimposed signal at is yu = a/pI^ia^i + \/Psht,2k s 2 + Tik, where 
Sj denotes the modulated symbols transmitted by Sj with the average power normalized, j = 1,2, and 
rifc is additive white Gaussian noise (AWGN) at R^, which is a zero mean complex-Gaussian RV with 
two-sided power spectral density of Nq/2 per dimension. During the second phase, i?& amplifies the 
received signal and forwards it back to the sources. Let be the signal generated by Rf,, then we have 
Xk = -JPrfikUk, where (3k is the amplification factor. In this paper, we analyze the variable-gain AF relay 
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IflBI . then /3fc = I p s /i^ife +p s h t) 2& + iVo ) is decided by the estimated instantaneous channel 
coefficients. 

The received signals by Si and S2 axe, similar due to the symmetry of the network topology, and thus, 
we take S\ as an example for analysis. The signal yi r received by S\ can be written as yi >r = ht 1 xk%k J r' n l' 
where m is AWGN at S\\ n\ and are i.i.d. RVs. According to (0]), h t ^ and ht 2 k can be rewritten 
as h t> ik = p e ht,ik + ^lfc an d h t 2k = Peht,2k + ^2fc> where difc and c?2fc are independent RVs due to the 
independence of h^ik and ht t 2k- Therefore, yi >r can be expanded as : 

Vl,r = VPrPsPkPlh t ,lkht,2kS2 (12) 
+ y/PrPsPk {Peh,lkd2k + Peh,2kd\k + d lk d 2 ^j S 2 (13) 
+ VPrPsf3kP 2 e h tjlk h t AkSi (14) 

+ VPrPsPk hpeht,lkdlk + ^lfe) Si (15) 

+ Vp^PkPeh t ,iknk + Vp^Pkdikn k + m , (16) 

where (fT2l) represents the useful information from S2 ; (IT3b represents the inter-interference from S2 caused 
by CEE; (fT4l and (ITSb represent the self-interference from Si itself which can be subtracted totally by self- 
canceling if CEE does not exist |9]. However, with CEE, Si can only reconstruct y/p r PsPkPlht,ikht,\k s i 
at the receiver. Thus, only ([141) can be subtracted totally, whereas the self-interference of (fT5T > is residual; 
(fT6l ) includes the amplified noise from i?& and the noise at S\. 

After self-canceling yJp T PshpV l t,ikht,ikSi from yi r , and then multiplied by h^ lk h* 2k to compensate 
the phase rotation, the processed signal yi at Si is : 

2/1 = K,ikK,2k (yi,r ~ VPrPsPkPlht,ikht,ikSi) ■ (17) 
The transmitted information S2 can be recovered by maximum likelihood detection: 



S2 = arg mm 

Sn&A 



Vl ~ VPrPsfiPe 





2 




2 


2 




ht,ik 




ht,2k 


4 


5 


(18) 



where ||-|| represents the Euclid-distance, A is the alphabet of modulation symbols, and §2 is the 
recovered signal. 

According to (fT71) . the instantaneous received SNR 71 at Si can be written as : 



IprlpsPl 



71 



h+ik 



h+2k 



(5ip^ s p 2 a 2 D + ^ r p2 + 



t,ifc 



+ {^r^sP 2 e Cr 2 D + Vs) 



t,2fc 



(19) 



where ^ s = p s /No, ip r = p r /No, p e is the CEE coefficient, and the CEE variance a 2 D = 1 — p e . 
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In high SNR regime, p e — > 1 and cr^ = (1 — p e ) — >• 0, then the item Sf/vV's !) + Vvcf) + 1 in the 
denominator of (fT9l ) approaches 1, which can also be ignored when SNR approaches infinity 0. 
Therefore, 7! in high SNR regime can be simplified into : 



7l 







2 _ 




2 


a 


ht,ik 


b 


ht,2k 





HAk 



+ b 



where 



VvPe 



ht,2k 

pi 



71 in (l20l) is greater than that in (fT9l ). whereas they match tightly in high SNR regime. Therefore, we 
use 71 in (f20l > for asymptotic analysis in the followings. 



hlp r 1p s p 2 e a 2 D + VrPe + 



(20) 



(21) 



S. Distribution Function of the Received SNR 



The distribution of 71 in (1201) is decided by % 



,ik 



ik 



and %2k 



t,2k 



, which are determined 



by 7s,ife = h Si \k and ^ S) 2k = ^s,2fc according to Lemma 1. Furthermore, the distribution of 7 Si ifc 
and 7 Sj 2fc can be obtained by the above selection criterion. After some manipulations, we have 
Theorem 1: With the definition that : 

„2^2 „ „ c„/, „/, „2^2 1 „/, J2 



a A£e_ 1 + IprPe^D & A Pe _ ^r^sPe^D + VVPe + *Ps 



a tprPe b 



Ipr^sPt 



the cumulative distribution function (CDF) of ^\ is : 

N-l N-l 



F 71 (z) 



m=0 ?i=0 



TV - 1\ /TV - 1\ (-l) m (-1) 



2m + 1 2ra + 1 



(/1 + /2 + h + h) 



where 



fi = 2v abz exp (— (a + 6) z) If 1 ( 2zV abj , 

2(n + I) a 



16n 2 ab 



x if 



2(n+l) [(2n + 1) (1 - pi) + 1] 
I 2z 



z exp 



(2n+l)(l-pi) + l 



+ b)z 



2 (n + 1) ab 



(2n + l)(l-pi) + iy 

/ 8m 2 ab 

y (m + 1) [(2m + 1) (1 - pa) + 1] 

x i^i ( 2z 



z exp 



a + 



2(m + l)6 



(2m + 1) (1 - p 2 ) + 1 



2ab(m + 1) 



(2m + l)(l-p 2 ) + l 



(22) 



(23) 



(24) 



(25) 



(26) 



s 



16m 2 n 2 ab 



y (m + 1) {n + 1) [(2n + 1) (1 - pi) + 1] [(2m + 1) (1 - p 2 ) + 1] ' 

/ / 2(n + l)a 2(m + l)6 

x ex P - 77^ 7777 7 + 



(27) 



(2n + l)(l-pi) + l (2m + l)(l-/9 2 ) + V 

^ (a I ab (rn + 1) (n + 1) \ 

x ifi 4z4 / — \ '— . (28) 

X \ V [(2m + l)(l-p 2 ) + l][(2n + l)(l- Pl ) + l)J 

And K\ (■) is the first order modified Bessel of the second kind 11231 . (JT) is the binomial coefficient, 
and pj, j = 1, 2 satifies ([8]) in Lemma 1: p,, = 1 if py = 1, and P j = p e Pf j if Pf < 1- 
Proof: The proof of Theorem 1 can be found in Appendix B. ■ 
Due to the symmetry, it can be proved similarly that the CDF of the received SNR 72 at S 2 have the 
same form as 71, and their PDFs can be obtained by differentiating the CDFs. 

C. Asymptotic Performance of Average Symbol Error Rate 

For many common modulation formats, the average SER can be obtained by ifTBl : 



SER = aE 



a I' ft 2 



F~[ — ] e~dt , (29) 





where 7 is the instantaneous received SNR, Q (•) is Gaussian Q-Function ll23l . and a = 1, B = 2 for 
BPSK, a = 1, B = 1 for QPSK, a = l/log 2 M, j3 = log 2 Msin 2 (vr/M) for MPSK (M > 4). 
Applying Theorem 1 and d29l , the exact average SER of S\ can be obtained by IT241 (6.621.3)]: 

Jo ^u[Px)dx {a + pr+v r(M + l/2) 2'^ 2'« + /37 

where V (■) is Gamma function, and F (•) is Confluent Hypergeometric function 11231 . However, the exact 
form is too complicated to analyze the performance, thus we resort to the high SNR analysis [22J. 

Theorem 2: The asymptotic performance of SER in high SNR regime can be obtained in two different 
cases according to whether the CSI is outdated or not. 

• When the CSI is not outdated, i.e., the delay coefficients satisfy p^ = pf 2 = 1 , and the CEE 
coefficient p e is arbitrary, the average SER of Si in high SNR regime is: 

a (2N)\ f fl + ip r p 2 e a 2 D \ N f 5^ sP 2 e a 2 D + VvPe + Vs X ' ' 



SEB * -w~nt\\ api ) + v — e i^r — ) ) ' (31) 

where a and 8 are decided by the modulation format in d29l : ip s = p s /Nq and ip r = p r /No\ N\ is 
the factorial of iV; a 2 D = 1 — p e . 
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When the CSI is outdated, i.e., pf x < 1 or pf 2 < 1, and p e is arbitrary, the average SER of Si in 
high SNR regime is: 

^ ( 1± ^#) N T, (-T ( N : l ) 7, p l „, „ (3 2) 



2/3 V ^rP\ ) ^ ' V n ) (2n+l)(l- ft ) + l 
2/3 V VvV^ y J^ l J V m y(2m + l)(l- 



(2m + l)(l-p 2 ) + l ' 



where pj, j = 1, 2 satifies ([8]) in Lemma 1: pj = 1 if pj = 1, and = p e p/j if P/^ < 1. 
Proof: The proof of Theorem 2 can be found in Appendix C. ■ 
Similarly, the average SER in high SNR regime of S2 can be obtained from (fJU and (1321 ) by permuting 
a with 6. 

D. Performance Analysis of Diversity Order and Coding Gain 

Diversity order d = — lim (log SER^ / log V't) ll22l . where ip t = (2p s + p r ) /N = 2i\) s + ip r , is 
an useful metric to describe the asymptotic performance of SER, i.e., greater diversity order means the 
curve of SER attenuates more quickly. 

Theorem 3: According to the definition of diversity order, the diversity order is : 



d 



N, if the CSI is not outdated, i.e., pj 1 = pf 2 = 1; 



(33) 



1, if the CSI is outdated, i.e., pj x < 1 or pf 2 < 1. 
Proof: Assuming ip t = (2p s +p r ) /N = 2tp s + tp r , tp s = p s /N = Xip t , and ip r = p r /N = 
(1 — 2 A) ipt, the diversity order can be obtained by Theorem 2, and the fact that p e — > 1 and a 2 D = 
1 — p e — > when SNR approaches infinity. ■ 
Theorem 3 reveals that the diversity order is N if and only if the CSI is not outdated. Once the CSI 
is outdated, i.e., the delay exists, the diversity order reduces to 1, whereas CEE has no impact on the 
performance loss of diversity order. 

However, both delay effect and CEE can reduce the coding gain, which is the shift of SER curve, e.g., 
different delay coefficients pj. and CEE coefficients p e will result in different pj is Theorem 2, and thus 
the coding gain is different. 



IV. Simulation Results and Discussion 

In this section, the average SER of bidirectional relay selection with imperfect CSI is studied by Monte- 
Carlo simulations, and the analytical performance provided by Theorem 2 is verified by these simulation 
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results. Due to the symmetry of the network, the following results only concern about the average SER 
of Si. All the simulations are performed with BPSK modulation over the normalized Rayleigh fading 
channels. For simplicity, we assume that sources and relays have the same power, i.e., p s = p r = P , 
and the x-axis of the following figures is SNR = Pq/Nq in decibel. To better understand the impact of 
imperfect CSI, we discuss four different situation, i.e., perfect CSI, only delay effect, only CEE, and 
both delay effect and CEE. 

In Fig. 1, we compare the simulated and the analytical SER of bidirectional relay selection with 
perfect CSI for N relays, i.e., pf t = pj 2 = 1 and p e = 1. This figure shows that increasing the number 
of available relays can reduce the average SER, because the diversity order is N when the CSI is perfect. 
This figure also shows that the asymptotic analytical SER given by Theorem 2 is the lower bound of the 
simulated results due to the fact that 71 in d20l ) is greater than that in ( fT9l ), whereas both the analytical 
and the simulated results match tightly in high SNR regime. 

In Fig. 2, we analyze the impact of delay on the SER performance without CEE , i.e., p e = 1. For 
simplicity, we assume p* = pf 2 = pt and N = 4. The figure reveals that the diversity order degrades 
to 1 once pj < 1 regardless of N. Although the diversity order is 1 once pf < 1, yet the coding gain is 
different for different pj. Comparing the curves of pf = 0.9 and pf = 0, the coding gain gap between 
them is approximately 6dB in high SNR regime. Besides, the performance at moderate SNR is different 
for different pj, i.e., greater pf has better performance at moderate SNR. For example, at moderate SNR, 
i.e., range from 8dB to 16dB, the slope of the SER curve of pf = 0.9 is greater than 1, while the slope 
of pf = at the same range is 1. The performance at moderate SNR can be analyzed by the exact 
expression of SER and Maclaurian Series ll22l . 

In Fig. 3, we study the impact of CEE on the SER performance without delay, i.e., pf = 1 and 
p e = P / {P + Nq), where P is the power of the training symbols Ell . P can be greater than the power 
of the data symbols Pq to obtain better performance of channel estimation, thus we simulate the situation 
of P = Po,2Po,APo and 00 (P = 00 means no CEE), respectively. With CEE, the diversity order is 
invariant, which is the same as the number of relays. However, compared with the curve of P = 00, 
there exists coding gain loss caused by CEE, and the loss could be reduced by increasing the power of 
training symbols Pq. As Fig. 3 illustrated, the coding gain loss in high SNR regime is about 5dB when 
P = Po, but it reduces to 2dB when P = 4P . 

In Fig. 4, the joint effect of delay and channel estimation error is considered and compared with the 
cases of only delay effect, only CEE, and perfect CSI. The results also indicate that delay will result in 
the diversity order loss and the coding gain loss, and CEE will merely result in the coding gain loss. 
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With both delay and CEE existing, the SER performance is the worst, which matches tightly with the 
analytical result in high SNR regime. 

V. Conclusions 

In this paper, we analyzed the performance of bidirectional AF relay selection with imperfect CSI, i.e., 
delay effect and channel estimation error, and the asymptotic analytical expression of end-to-end SER 
was derived and verified by the computer simulation. Both analytical and simulated results indicate that 
delay effect results in the coding gain loss and the diversity order loss, and channel estimation error will 
merely cause the coding gain loss. 

Appendix A 
Proof of Lemma 1 

At the case of pf. = 1, we have htji = h s ji by £[]), and thus htji = h s ,jh which is a special case of 
f7]) when pj = 1 in Lemma 1. 

At the case of < 1, by (Q}, © and ©, we have : 



ht,ji = hji + e = Pfrhgji + J 1 - pj.Sj + e = pf 3 p e h s ji + pf.d + J 1 - pj.Ej + e , (34) 
where d, £j, and e are independent zero mean complex-Gaussian RVs with variance of a 2 D , a 2 ., and 



a e , respectively. Thus, pf.d + wl — pj Ej + e is a zero mean complex-Gaussian RV with variance 

of „2 _l_ „2 ^ ^2 , ^2 „,u^U „„„ k» ( 1 „2 „2^\ ^2 



Pf °D + y- ~ Pf J a e 3 + °e > which can be simplified into ^1 — pf.pij cr| by the relationship 
of variances (fS),©. Then, pf.d + yj\ — p 2 j £j + e can be written as yj\ — pj.pivj , where vj is an 
independent RV with zero mean and variance of cr? . Defining pj = p e Pf 3 . formula f7]) in Lemma 1 

2 

are 





2 




htji 


and 7 s jj = 


hs,ji 



is proved. Thus, htji and h s ji are jointly complex-Gaussian, and jt,ji 
correlated exponential distributions, then the joint PDF f^ t ^ (y, x) is given by (25 

fj tj i,%ii (Vi x ) = 7 -\ ex P f ~ 7 - ' \ I 7 o I 7 V \ I • ( 35 ) 




3 ) h 

And now, the conditional probability of (O in Lemma 1 can be proved by |[26l : 



h, 3 ,\%,, (y\ x ) = j (a) — > ( 36 ) 



where f % ji (x) = exp (-ff/ffj) l°\ ■ 
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Appendix B 
Proof of Theorem 1 

A. distribution of *{ s ik and ^ s ^k 

Following the similar steps of lfl4l , the CDF of 7 Sj ifc can be expressed as : 

F % lk (x) ( = } N Pr {j Sili < x,k = i} (37) 

X 

= N / fj BU (y)Pr{% :li < % : 2i\%,u = y}Pv{k = < J s , 2 i,ls,li = y} dy 



+ N J f% Ai (y) Pr{7s,ii > 7 s ,2i|7s,ii = y) Pr i k = > %,2i,%,ii = y} dy 



where (a) in d3"7T ) is satisfied due to the symmetry among the N end-to-end paths, and (b) is satisfied 
by dividing the union event into two disjoint events, i.e., 7 Sj ij > j St 2i and 7 Sj ij < 7 Sj 2i- According 
to the selection criterion (fTTT > and order statistics of independent RVs |[27l : Pr {min (xl, x2) < = 
1 - (1 - F Xl (z)) (1 - F X2 (z)), and the fact that F % li (z) = F % _ 2 . (z) = 1 - exp (-z/o|), we have 

Pr{£; = i|7 SjU < 7 Si2 i,7s,ii = y} = Pr {min (7^, 7 S)2i ) < y} = ^1 - exp ^-^|^ . (38) 

Similarly, the conditional probability Pr{fc = i|7a ii > 7s,2i>7s,ii = 2/} can be achieved. Therefore, sub- 
stituting (1381 ) into (l37l ). lfe (x) can be written as : 

N-l 



— exp ( — ^ ] ( 1 - exp ( — | ) ) ./: j dy (39) 



N-l 

oo 



+ N [ exp ( — ^ I ( / ^2 exp ( — ^ J J ( 1 - exp ( 1 I I dy . 

Applying binomial expansion (1 - xf = £f =0 (?) and ^ En=o (V) ("^V ( n + X ) = 1 

in |[24l (0.155.1)], F* lk (x) can be rewritten as : 



•n—n \ / 



n J In + 1 



(x\ n / 2 (n + 1) x 
2 M TT ex P 2 



(40) 



n=0 

where cr? = p" 1 , and it can be proved similarly that the CDF of 7 Sj ifc have the same form, and their 
PDFs can be obtained by differentiating the CDFs. 
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B. distribution of 71 



At the case of pf 1 < 1 and pj 9 < 1, and by Lemma 1 and J exp (-ax) I (Py/x)dx = (1/a) exp (/3 2 / (4a)) 



in 121 (6.614.3)], we have : 



2V-1 



n=0 



JV— 1 



1 



n J 2n + 1 



exp ( —xjo\ 



o~7 



(41) 



+- 



2n/a? 



2(n + l)a?/<7? 



(2n + 1) (1 - pf) + 1 6XP ^ (2n + 1) (l - pf) + l y 
The CDF of 7^1^ can be obtained by integrating the PDF, and the distribution of jt,2k can be obtained by 
substituting p\ with p^. Let Qi and Q2 represent d^t.ik, and 67^2/c respectively, and the distribution of Q\ 
and ^2 can be obtained by fy (y) = (1/m) fx (x/m), and Fy (y) = Fx (x/m) when Y = mX (m > 0) 
1126*1 . Thus, the CDF of 71 can be written as : 



F 7l (z) = Pr <i — < z 



fii + ^2 

Pr {(0 2 - z) Oi < z^ 2 , ft 2 > z} + Pr {(^ 2 - 2) fii < ^2, ^2 < z] 



(42) 



1 - y /n a (x + z) 





1 - \z + 



1-^ 



.c.r 



x — z 



fu 2 (x) dx 



Substituting J exp (-mx — nx x )<ix = 2\J~njmK\ (2y/mn) in [24 , (3.324)] into (1421 ). Theorem 2 can 



be proved when using NJ2n=o CVK -1 )"/ (n + 1) = 1 in (23 (0.155.1)]. 

At the case of p/ x = 1 or pf 2 = 1, it can be proved in a similar way that the CDF of 71 can also be 
expressed as the formula (f23l in Theorem 1. 



Appendix C 
Proof of Theorem 2 

In high SNR regime (ip r ,ips —> 00), p e — > 1 and a, 6 — >• 0. By applying the Bessel function 
approximation for small x — > 0, K 1 (x) « 1/a; [23] in Theorem 1, we have : 

'JV-l\/jV-l\ 1 1 
n 



N-1N-1 



f 71 (*) « 1 - iv 2 E (-*: 



rn+n 



m=0 n=0 



111 



n 



exp (— (a + b) z) -\ exp 

n + 1 



2m + 1 2re + 1 
2(n+l)a 



(43) 



( 2n + 1) (1 - px ) + 1 



+ 6 U 



m 

+ — -exp 

m + 1 



a + 



2(m + 1)6 



;2m + 1) (1 - p 2 ) + 1 
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m n 

H exp 

m + ln+1 



2(n + I) a 



+ 



2(m + l)6 



(2n + 1) (1 - pi) + 1 (2m + 1) (1 - /9 2 ) + 1 



At the case of pf ± = pj 2 = I, F 7l (z) can be rewritten as : 

AT-l 



F 7l M 



>-*£<-i>"( „ j 

n=0 v 7 



n / 2n + 1 



exp (— az) H -j— j- exp (—2 (n + 1) az) 



(44) 



AT-l 
m=0 

Furthermore, we have : 

AT-l 



N - 1\ 1 



m 



/ 2m + 



exp (— bz) H exp (—2 (m + 1) 6z) 

m + 1 



iv£(-iy 



n=0 
AT-l 



AT - l\ 1 



n / 2n + 1 



exp (—az) H exp (—2 (n + 1) az) 

n + 1 



(45) 



/AT_1\ 1 
^E^H exp(-2(n+l)a.) 
/ — ' \ re / n + 1 



n=0 
N—l 



/AT — A 1 

+ iV ^ (-1)™ — - - [exp (-az) - exp (-2 (n + 1) az) 
n=o \ n y n + 



(a) 



(6) 



1 — [1 — exp (—2az)] N — exp (—az 



N-l\ 1 



n 7 2n + 1 



E 

p=i 



(-2az) % 
p\ 



N 



£ 

p=0 



AT-l , 

iv^(-ir( 

n=0 ^ 

n=0 ^ 7 p=l 



(1 - exp(- (2n + l)az)) 



(-az) J 



® 1 - i (2az) w + ((az)") , 



where (a) in © is achieved by the fact that ( A/ " 1 ) ^ = ( n ^) and (1 - x) N = J2k=o ( ^ ) (-1) (b) is 
achieved by Maclaurian Series of exp (x) = Y^^Lo [x p ) / (pi); (c) is achieved by Y^£L (~^-) k Cn)^" 1 = 
0, (1 < n < N) in (21 (0.154.3)]. Therefore, 



F 71 (z) « [5 (2a)" + i (26) V* 



(46) 



Finally, @D in Theorem 2 is proved by @2), ([29]), @, and / f 27V exp (-i 2 / 2 ) eft = 2( Ar " 1 / 2 )r (1/2 + N) 

o 

y/n/2(2N)\/(2 N N\) in JM (3.326.2)], where T (•) is Gamma function l23l 
At the case of pt t < 1 or pj 2 < 1, by exp (x) « 1 + x, we similarly have : 



N ^ //V-1\ 



2 -Pi 



n=0 
AT-l 



n 7 (2n + 1) (1 - pi) + 1 



(47) 



-6ziV £ (-1)' 



m=0 



JV-1 



2-P2 



m / (2m + 1) (1 - p 2 ) + 1 



Then, (|32l ) in Theorem 2 is proved by (1221) and d29l ). 
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Fig. 1. Analytical and simulated SER with perfect CSI, with different N, pf 1 — pf 2 =1 and p e = 1. 
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Fig. 2. Analytical and simulated SER of Si with delay effect, with different pf, and p e = 1, N = 4. 
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Fig. 3. 



Analytical and simulated SER of Si with estimation error, with different p e and pf ± = pf 2 — 1, iV = 4. 
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